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Introduced as a model for hyperchaos, the generaUzed Rossler system of dimension A*' is obtained 

by linearly coupling TV ~ 3 additional degrees of freedom to the original Rossler equation. Under 

variation of a single control parameter, it is able to exhibit the chaotic hierarchy ranging from fixed 

points via limit cycles and tori to chaotic and, finally, hyperchaotic attractors. By the help of a mode 

transformation, we reveal a structural symmetry of the generaUzed Rossler system. The latter will 

allow us to interpret number, shape, and location in phase space of the observed coexisting attractors 

Qs ' within a common scheme for arbitrary odd dimension N. The appearance of hyperchaos is explained 

^^ I in terms of interacting coexisting attractors. In a second part, we investigate the Lyapunov spectra 

^\ . and related properties of the generalized Rossler system as a function of the dimension A^. We 

T-H ' find scaling properties which are not similar to those found in homogeneous, spatially extended 

f— I , systems, indicating that the high-dimensional chaotic dynamics of the generalized Rossler system 

I— { ■ fundamentally differs from spatio-temporal chaos. If the time scale is chosen properly, though, a 

i—i ' universal scaling function of the Lyapunov exponents is found, which is related to the real part of 

W-, , the eigenvalues of an unstable fixed point. 
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Chaotic dynamics has been intensively investigated with the help of simple low-dimensional models like the Lorenz 
or the Rossler system [0 . Because of the restricted phase space, only low-dimensional chaotic motion is observed in 
^\ , these systems. Low-dimensional chaos can also be observed in nature 0, although the underlying dynamical systems 
have an infinite number of degrees of freedom. One may wonder: What is the connection of this low-dimensional 
chaos to the truly high-dimensional states, that may arise out of chaos under variation of one control parameter? How 
i-pj '■ does the transition from the low-dimensional to the high-dimensional states take place? 
I ' Consider the attractors of dissipative ordinary differential equations under variation of the dimension N of phase 

2 ■ space. Rossler ||^ has postulated a chaotic hierarchy where more and more qualitatively new forms of complex motion 
develop with increasing dimension: In one dimension, only stable fixed points are encountered. In two dimensions, 
^ also periodic orbits can exist. In three dimensions, the possibility of quasiperiodicity and chaos arises. Thus, at the 
lower end of the hierarchy, the well-known low-dimensional dynamical states of motion can be found. Somewhere high 
up in the hierarchy, things like turbulence or noise may be located. One crucial question is, how one can distinguish 
and characterize the higher steps of the hierarchy. A provisional classification can be made in terms of Lyapunov 
5^ ] exponents [p[. For dynamical states with more than one positive Lyapunov exponent, the term hyperchaos has been 
coined 0. In order to study the full chaotic hierarchy, Baier and Sahle M have introduced a class of model equations, 
the Generalized Rossler System (furtheron called GRS). Starting from the Rossler system as one of the simplest and 
best understood nonlinear ordinary differential equations that exhibit chaos, the GRS is obtained by linearly coupling 
additional degrees of freedom to the original Rossler system. The structure of the GRS of dimension N is that of an 
TV — 1 dimensional linear subsystem that is coupled to one nonlinear variable. As will be shown, the GRS preserves 
essential characteristics of the Rossler system while extending it to a phase space of arbitrary dimension N. This 
allows to study the influence of the discrete parameter N, without having to compare completely different systems. 

Baier and Sahle W demonstrated that the GRS does, indeed, show hyperchaotic dynamics with an increasing 
number of positive Lyapunov exponents for increasing N. The GRS realizes one possible path through the complete 
chaotic hierarchy from a stable fixed point via periodic orbits and chaos up to hyperchaos. We ||^ have introduced a 
mode transformation of the GRS based on the numerical solution of the linear subsystem. In the present paper, we will 
restate the mode transformation based on a semianalytical solution of the linear subsystem. The mode transformation 
will then be used to analyze the dynamics of the GRS with arbitrary dimension N in phase space. In Section O, 
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the GRS is introduced and general properties of the GRS are discussed. In Section III, we deal with the mode 
transformation and the concept of structural symmetry, that later on will allow us to understand the number, the 
form, and the location of the coexisting attractors of the GRS. Subsequently, we numerically investigate the dynamics 



of the GRS in phase space for the cases N = 5 and A^ = 7 in Section IV . We interpret the observed dynamics within a 



general scheme that enables to predict the structure of the attractors for higher N. In Section ^ then, we investigate 
the GRS in the hyperchaotic state for different values of the dimension N. We report scaling properties of the number 
of positive Lyapunov exponents, the Lyapunov dimension, the metric entropy, and the Lyapunov spectra as a function 
of N. The limit N — > cx) is discussed in view of the literature on that subject. 



II. GENERAL PROPERTIES OF THE GENERALIZED ROSSLER SYSTEM 



The GRS is given by 



— x(t) = Ax{t) - a;Ar(i)ejv-i 
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It is composed of an A'^— 1 dimensional linear subsystem x and one nonlinear variable xn- Qualitatively, the mechanism 
of instability of the GRS is the same as that of the Rossler system for all dimensions N. The positive feedback or 
autocatalytic process that is controlled by the parameter a causes an expansive dynamics of the GRS around the 
origin. As long as xn-i remains well below the threshold d of xn, x^ adiabatically follows its equilibrium value 
e/b{d — xjv-i) and does not influence the linear subsystem x appreciably. When xn-i comes close to or exceeds d, 
xjv will start to grow rapidly, thereby folding the system back to a state of lower amplitude (via ejv_i in Eq. (^). 
This time development of xn of long intervals of small amplitude interrupted by short spikes leads us to call xn the 
nonlinear trigger. 

One of the three parameters e, b, and d can be eliminated by rescaling the amplitude of (x, xat). In the form the 
equations are given in dl^) , (0) , they are scaled in time in such a way that the (angular) eigenfrequencies of the linear 
subsystem lie in the interval [0,2] for all values of N (see Section HI). 

The divergence of the GRS is a + b{xN-i — d), independent of N. Thus, the GRS is dissipative, if (xjv-i) < d — |- 
(the angle brackets denote the time average). This is true for all parameter values to be considered in the following. 

For odd N, the linear subsystem x can be transf orm ed into ~ harmonic oscillators that are coupled only via the 
nonlinear trigger xn, as will be shown in Section III (compare also [p|). For even N, the linear subsystem can be 
transformed into ~^ oscillators together with one variable, which simply grows exponentially. That means, there is 
one positive real eigenvalue in addition to -^^2^ pairs of com plex c onjugate eigenvalues. This leads to a qualitatively 
different dynamics. See the remark at the end of subsection IIIB for clarifying this point. In the present paper, we 
restrict ourselves to the case of odd N. 



III. THE MODE PICTURE 



In the following section, we develop the mode picture of the GRS. First, we will present the solution of the linear 
subsystem. Herewith, we can transform the GRS into a mode picture, where the linear subsystem x consists of 
harmonic oscillators, which are coupled only via the nonlinear trigger variable x^. In the mode picture, we are in the 
position to explain the idea of the structural symmetry, which will be of central importance for understanding the 
dynamics of the GRS in subsequent sections. Finally, we show the fixed points together with their stability properties. 



A. Semianalytical solution of the linear subsystem 

Consider the linear subsystem x (Eq. (|l|)). Here, we restate the equations in component form as a hnear chain with 
boundary conditions: 

c„+i, ne {!,... ,N ~1}, 

(3) 

To solve the linear subsystem x, an exponential ansatz is used: 
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xan+i = cos((2n + l)fc + ^)e*"('=)* 

, forne {0,...,^} 
X2n = csin(2nfc + (y5)e*^('^^* 

, forne {0,...,i^} 



(4) 



where k is the wave number of the eigenmodes, fl{k) the corresponding frequency with a specific dispersion relation. 
Substituting ansatz (^) into chain (||) yields 



c = —I 



n±„i = 2 sin(fc±,„) , TO e {1, . . . , \^} , 

7, I 27n — l __ 1 



(5) 



These (A^ — 1) eigenmodes are a complete solution of the {N — 1) dimensional linear subsystem x for odd iV. The 
solution is analytical up to a complex correction Lp±„i to the allowed wave numbers k±jn, which is determined by a 
transcendent equation. For a = 0, this equation can be solved giving 

k = ^f^TT , where to e {1, ... , ^} . ^^^ 

For a < 1, there are Z^ pairs of complex conjugate eigenvalues (ri,„, il_m), where ilm = —fl-m- The corresponding 
pairs of eigenmodes represent harmonic oscillators with angular frequencies Um = Re fim and autocatalytic coefficients 
a,„ — —2 * Imi7ni. The complex corrections 4'±m to the values of /c„ for a ^ can easily be expanded about a = 0. 
The expansion is up to the order of a^ 



^±m = MC0s(22^7r) 

+ ^«'sin(±2l^7r) + 0(a3). 



(7) 



In the considered range of a G [0,0.35], this is an excellent approximation (however, in the subsequent analysis, we 
have numerically calculated the eigenvalues and eigenvectors of A in order to minimize errors). The imaginary part 
of ifm determines the autocatalytic coefficient am- It is linear in a up to O(a^). The real part leads to a correction of 
LUm- But for a £ [0, 0.35], this correction is very small, i.e., the frequencies of the oscillators are almost independent of 
a in the considered range. From equation M), we infer that the frequencies of the oscillators are bounded (0 < a; < 2 
for all N) and we expect the autocorrelation time t^ in the case of an unstable dynamics to be independent of N. 
Additionally, the velocity of signals of frequency u traversing the linear subsystem equals 2 cos k. Therefore, the time 
Ts a signal takes to traverse the linear subsystem (from xi to a;(7v-i) or back, since the linear subsystem allows for a 
bidirectional transport of signals) is expected to be approximately N/2. 

B. The mode transformation 

We have solved the linear subsystem x. The next step is to transform it into its eigensystem y. The transformed 
equations for y are completely decoupled (as the matrix A is transformed into a diagonal matrix). Thus, we can treat 
the different oscillators independently. The idea of the following steps of the transformation is to bring the oscillators 
that correspond to pairs of complex conjugate eigenvalues into a form that is as similar to the oscillator of the original 
Rossler system as possible. The oscillators are of the form 

dt\ y2rn{t) y ~ V -^f}„^ ) \ y^^it) 



It can be considered as the principal value decomposition of the following Rossler-like oscillator: 
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In the following, these oscillators are called oscillator {z2m-i, ^2m)- For each variable of the linear subsystem x, one 
can choose one complex factor of normalization. It can be utilized to make the entire transformation real (via phase 
factors) and to give the coupling to the nonlinearity xn - which is not touched by the transformation - a certain form 
(via amplitude factors) . We choose this coupling in such a way that the resulting equations are 



where 
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as the divergence of the GRS is not affected by the transformation. 



c denotes the vector of coupling constants of xn to the components of z. If the complete transformation is 



= V-'bn- 
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x), then the vector of coupling constants c is c 

^ 0, (cosfc w-i )^ ) . This mode transformation converts the GRS into a system of 



called U (i.e. z = U^^ 

c — -^1 0, (cos/ci)^, 0, (cosfc2 

oscillators with frequencies ujm and autocatalytic coefficients am, that are coupled solely via the nonlinear trigger 
xn- The coupling has a special form: the oscillators couple to the trigger only with the sum of their variables Z2rm 
i.e., they couple to xn identically. The coupling of xn back to z is different for the different oscillators. In the 
following, we call the original form of the GRS the Baier-Sahle picture and the transformed form the mode picture. 
In Fig. |l|, the structures of the Baier-Sahle picture and the mode picture are compared. Each linear degree of freedom 
is represented by a circle, each nonlinear degree of freedom by a square. Couplings are shown as connecting lines. In 
the Baier-Sahle picture, the similarity of the GRS to a spatially extended system is most obvious; the autocatalytic 
process and the nonlinear trigger xn can be considered as the boundary conditions of a homogeneous linear chain. 
On the other hand, in the mode picture, the dynamics of the GRS can be understood as the interaction of different 
oscillators, that are coupled only via one nonlinear trigger. As the trigger xn influences the dynamics of the linear 
subsystem only during the presence of the spikes, the effect of coupling is restricted to these short intervals in time. 
This allows to obs erv e the dynamics of the different oscillators independently in the respective projections of phase 
space (see Section IV). 

For even values of N, one finds one positive real eigenvalue of the linear subsystem for a > 0. In the mode picture, 
this corresponds to one exponentially growing degree of freedom, which is coupled to the nonlinear trigger xn in the 
same way as the even coordinates of the oscillators. The contribution of the nonlinear trigger to the time derivative 
of this exponential mode is always negative. Thus, if once the exponential mode is pushed to a negative value by the 
nonlinear trigger it will escape to — cx). This mechanism, in general, leads to a global instability of the GRS for even 
N. 



C. The structural symmetry 



The mode picture reveals the GRS to consist of harmonic oscillators (z2m-i,Z2m), that are identically coupled to 
the nonlinear trigger xn (compare Fig. |^). A perfectly symmetric system of the same structure as that in Fig. |l|(b), 
i.e., one with identical coefficients for each oscillator, would be symmetric with respect to any permutation of the 
oscillators. The invariant manifolds of such a system, specifically, the orbits and attractors, would have to obey this 
symmetry. For the invariant manifolds, there are two possibilities. Either they show the full symmetry in themselves, 
i.e., they are symmetric with respect to any permutation of the oscillators. Or they have coexisting mirror images. 



such that the union of them obeys the full symmetry. In the GRS, the symmetry is broken merely by the difference in 
the frequencies Wm, autocatalytic coefficients am-, and coupling constants c„. Nevertheless, the symmetry is preserved 
as a qualitative feature of the dynamics, as will be shown later on. We call this property structural symmetry. For 
arbitrary odd A'^, the structural symmetry will be utilized to understand and predict number, shape, and location of 
the attractors of the GRS in phase space for moderate values of the autocatalytic coefficient a. The first manifestation 
of the structural symmetry can be seen in the stability properties of the fixed points of the GRS. 

D. The fixed points 

The fixed points of a dynamical system are the pivots around which the system evolves. Thus, it is essential to 
investigate the stability properties of the fixed points, if one wants to develop any understanding of a dynamical 
system. For odd N, the GRS has two fixed points which are the same in all odd dimensions N, in the sense that the 
common linear coordinates and the nonlinear coordinate are identical for any two GRSs of different N. The two fixed 
points are (for odd N) 

(1,2) _ (1,2) _ _ (1,2) 

X 1 Xq —— . . . JJ A 
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Here, the minus signs belong to the index (1). The fixed point {x'-^\xj^ ) lies close to the origin of the system (for 
e = 0, it would be the origin). Accordingly, the stability analysis of this fixed point yields, in good approximation, 
the eigenmodes of the linear subsystem. The additional A^-th eigenvalue is strongly attractive. It corresponds to 

the exponential decay of xjv to , ,_'^ and its value is approximately —bd, as is to be expected. On the unstable 

manifold of this fixed point, the GRS expands until it is folded back by a trigger event of xn- The presence of the 
nonlinear trigger xn slightly stabilizes the eigenmodes, such that they do no longer become unstable at a = 0, but 
approximately at a = 0.006. Interestingly, the oscillators that become unstable first are the ones that have a larger 
autocatalytic coefficient q:„ for higer values of a (see Fig. 0). 

At the fixed point (xf^^^xj^ ) , one can already perceive the fingerprint of the structural symmetry: for each 
oscillator of the GRS, there is an oscillatory instability of the inner fixed point {x^^\xj^ ). The symmetry breaking 



is responsible for the difference in the eigenvalues of the fixed point. It will be shown in Section IV that each of these 
instabilities of the fixed point {x.^^\x}f ) gives rise to a coexisting attractor of the GRS. The fixed point (x'^^^xjy ) 
governs the folding process of the GRS. However, none of the attractors ever comes close to this fixed point. 

IV. DYNAMICS OF THE GENERALIZED ROSSLER SYSTEM IN PHASE SPACE 

In the following section, the mode picture and the concept of structural symmetry will be utilized to develop an 
understanding of the dynamics of the GRS in the A^-dimensional phase space. At first, we discuss the case N = 5 
in some detail. It will be demonstrated that for values of a < 0.09 two attractors coexist. For higher values of a, 
they interact in several crises. Finally, the two attractors merge into one hyperchaotic attractor. This scenario of 
three parameter regimes of coexisting attractors, interacting attractors, and one large hyperchaotic attractor will 
then be shown to hold in the case A^ = 7 as well. Finally, we postulate that the dynamics of the GRS in arbitrary 
odd dimension N can be explained by extrapolating this scheme. For all numerical calculations, we have restricted 
ourselves to the variation of N and a. The other parameters remain constant at the values e = 0.1, b = 4, d = 2. 

A. Numerical study of tlie case N = 5 

In the mode picture, the five-dimensional GRS consists of two oscillators that are coupled via the nonlinear trigger. 
First, we present a bifurcation diagram (Fig. pi) and the corresponding Lyapunov spectrum (Fig. 0) under variation of 
the control parameter a (see Appendix for the numerical methods that were used) . The bifurcation diagram consists 
of two separate parts. They correspond to two coexisting attractors and have been obtained using different initial 
conditions. 



As one increases a from a = 0, the fixed point near the origin becomes unstable at about a — 0.0006. Shortly one 
after the other, the two attractors emerge in two subsequent Hopf bifurcations as limit cycles. These two periodic 
orbits constitute the basis for all further development of the attractors. Although they become unstable at some 
point, they exist up to the highest values of a. In the following, they will be referred to as the period-1 limit cycles 
of the two attractors. Attractor 1 develops as a result of the instability that corresponds to oscillator (zi, Z2) in the 
mode picture. Accordingly, attractor 2 develops out of (2:3, Z4). In Fig. g, two phase space projections of the period-1 
limit cycle of attractor 1 are shown: one onto oscillator (21,2:2) together with the nonlinear trigger variable X5 and 
the other onto oscillator (23, 24) togeter with 0:5. The amplitude of oscillator (23, 24) is negligible compared with that 
of oscillator (21, 22). Here, the dynamics of the GRS is completely dominated by (21, 22) together with X5. In general, 
the frequencies of the two oscillators are incommensurate. The two oscillators drift out of phase in each revolution. 
However, the spikes in x^ cause phase shifts in (23, 24) that resynchronize the two oscillators. This leads to a 1:2 mode 
locking. For higher a (i.e., above a = 0.05), oscillator (21,22) can no longer force the phase locking and the periodic 
orbit breaks up into a quasiperiodic one. As one increases a further, attractor 1 shows several periodic windows and, 
finally, low-dimensional chaos arises. In Fig. ^, the phase space projections of this low-dimensional chaotic form of 
attractor 1 are depicted. Oscillator (23,24) has still a much smaller amplitude than oscillator (21,22) (though its 
influence on the nonlinear trigger xq is no longer negligible). And, remarkably, the dynamics in the three-dimensional 
projection onto (21,22) together with 2:5 shows a close resemblance to the Rossler system. It, therefore, can be 
considered as a 'perturbed' Rossler system. 

The evolution of attractor 2 under variation of a is completely analogous, except that now oscillator (23, 24) plays 
the dominant role and (21, 22) is very small. In Fig.^, the period-1 hmit cycle of attractor 2 is depicted. It corresponds 
to the periodic orbit of attractor 1 shown in Fig. |^ Here, oscillator (23,24) is dominant, whereas oscillator (21,22) 
has a negligible amplitude. The phase locking between the oscillators is accomplished by phase shifts in oscillator 
(21,22) due to the spikes in x^ which are in turn triggered by (23,24). Again, the dynamics in the projection, here 
onto (23,24) together with 25, resembles that of the R5ssler system, which can best be seen on the low-dimensional 
chaotic orbit in Fig. g (compare the chaotic orbit on attractor 1 in Fig. g). 

The existence of the two attractors is a consequence of the structural symmetry: in a perfectly symmetric system 
(of the same structure as the GRS in the mode picture), each attractor, that is not in itself symmetric with respect 
to the exchange of any two oscillators (i.e., the transformation 21 <-*■ 23, 22 ^^ 24 in the case TV = 5), must necessarily 
have one or more mirror images, such that their union fulfills the symmetry. An attractor, where oscillator (21,22) 
dominates oscillator (23, 24), is obviously not symmetric with respect to the above transformation. Therefore, in the 
perfectly symmetric system, there would have to be another coexisting attractor, where oscillator (23,24) dominates 
oscillator (21,22). In the GRS, the symmetry is broken by the difference in the frequencies, autocatalytic constants, 
and coupling constants. However, as has been shown above, these two attractors do still exist. Moreover, their shape 
and location in phase space does, indeed, reflect the structural symmetry. Attractor 1 lies close to the subspace 
defined by 23 = and 24 = 0. Attractor 2 is situated close to the subspace defined by 21 = and 22 = 0. On attractor 
1, oscillator (21, 22) is dominant and, on the periodic orbits, the phase locking is kept up by phase shifts of oscillator 
(23,24). On attractor 2, accordingly, oscillator (23,24) is dominant and, on the periodic orbits, the phase locking is 
kept up by phase shifts of oscillator (21, 22). 

Up to a ci 0.09, the two attractors develop independently. They are located clearly separate in phase space. Then, 
a boundary crisis occurs: Attractor 1 collides with the separatrix to the basin of attractor 2 and ceases to exist as a 
stable dynamical state. It still attracts the trajectories out of its basin, but eventually every trajectory ends up in 
attractor 2. At a ~ 0.115, another crisis can be observed. In the bifurcation diagram (Fig. H), the dynamics on the 
resulting attractor resembles that of attractor 1. However, it is not at all obvious, whether it is, indeed, attractor 1 
or, maybe, the two attractors merged. In order to answer this question, we select an initial condition on the period-1 
orbit of "attractor 2" and observe how the GRS develops from there. For this purpose, we utilize the time-delayed 
feedback control method introduced by Pyragas M (see Appendix, for the utilized numerical algorithm). With the 
aid of the control, we force the GRS onto the period-1 orbit of "attractor 2" (Fig. 0) at a = 0.12 (see Fig. |^). It 
should be remarked here that, in a strict sense, one can no longer talk of the "attractor" after it has undergone a 
boundary crisis. However, as it still attracts the trajectories of its basin and the transient motion along them may 
be very long, we will, nevertheless, still call it attractor. The fact that this is possible (with vanishing control signal) 
demonstrates, that at least the period-1 limit cycle of attractor 2 docs still exist. Subsequently, we switched off the 
control. First, the trajectory left the unstable periodic orbit according to the largest Lyapunov exponent of that orbit. 
Then, it remained on attractor 2 for approximately 1000 revolutions. But, eventually, the trajectory changes very 
quickly to an attractor-1-type shape and, within the patience of the authors (for a relatively long time interval of 
At = 300000, i.e., about 30000 revolutions), it never again showed attractor-2-type behavior (Fig. noh. From this, we 
conclude that the attractor at a = 0.12 is, indeed, attractor 1 and attractor 2 has lost its stability. In order to arrive 
at this situation, two things must have happened: attractor 2 must have undergone a boundary crisis into the basin of 
attractor 1 and attractor 1 must have separated from the separatrix in an inverse boundary crisis. It seems virtually 



impossible that these two incidents occur at the same value of a. Thus, depending on the order they occur, there 
should be a small interval of either coexisting or merged attractors in between. There is some numerical evidence 
that, indeed, the second situation is valid. 

The development of the two attractors as a function of a can also be observed in the Lyapunov exponents. Figure ^ 
shows the Lyapunov exponents of attractor 1 except for the interval a £ [0.9, 0.115], where attractor 1 does not exist. 
The sequence of periodicity, quasipcriodicity with periodic windows, and chaos can clearly be seen. For a > 0.115, the 



GRS is hyperchaotic ||6|. In Fig. 11, the phase space projections of the hyperchaotic attractor at a = 0.3 are shown. 
Interestingly, the second Lyapunov exponent becomes positive in the vicinity of the second and the third crisis. From 
that point, the three largest Lyapunov exponents grow more or less linearly with a (exept for a small periodic window 
at a ~ 0.23). The GRS in this range of a can no longer be regarded as a perturbed Rossler system in any projection. 
Oscillator (zi, Z2) gradually ceases to be dominant and the dynamics of the GRS seems more and more irregular. At 
a ~ 0.17, the third Lyapunov exponent becomes positive. The Lyapunov dimension of the attractor is larger than 4. 
However, no qualitative change of the dynamics could be observed at this point, cither in the bifurcation diagrams or 
in the phase space projections. 

B. Numerical study of the case TV = 7 

We now proceed to show that the GRS does behave in a similar manner in the case N = 7. In the mode picture, 
the seven-dimensional GRS consists of three oscillators coupled via the nonlinear trigger 2:7. 

If one compares the bifurcation diagram (Fig. |l^) and the corresponding Lyapunov spectrum (Fig. |l3| ) with those 
of A'" = 5 (Figs. and H, respectively), one perceives a striking resemblance. There are three coexisting attractors for 
lower a {a £ [0, 0.05]), corresponding to the three oscillators of the seven-dimensional GRS. Two of these attractors 
vanish via boundary crises in an intermediate regime of interacting attractors (a G [0.06,0.115]). Finally, only one 
attractor remains, which is hyperchaotic with an increasing number of positive Lyapunov exponents. 

For a e [0, 0.05], each of the three attractors can be characterized by the dominance of one of the three oscillators. 
Each attractor develops as a result of the instability of the fixed point (x^-'^^xj^ ) corresponding to the respective 
dominant oscillator. This oscillator triggers 2:7, while the other two oscillators have a small amplitude. On attractor 
1, oscillator (21,^2) is dominant, on attractor 2, oscillator (23,24), and on attractor 3, oscillator (25,25). As long as 
the attractors are periodic (e.g., for a G [0.006, 0.04] on attractor 1) the mode locking is kept up by phase shifts in the 
two oscillators with the small amplitude as in the case N — 5. Most remarkably, on all three attractors, the dynamics 
is qualitatively equal to that of the Rossler system in the projection onto the respective dominant oscillator together 
with the nonlinear variable x-/, i.e., in a three-dimensional projection of the seven-dimensional phase space (compare 
Fig. ^ (a) and (b) for the period-1 orbit on attractor 2, and Fig. ^ (c) and (d) for a chaotic orbit on attractor 2). 
Again, the existence, form, and location in phase space of the three attractors can be interpreted as a manifestation of 
the structural symmetry. They appear as mirror images of each other under the symmetry operation of an exchange 
of the dominant oscillator with one of the other oscillators. For example, attractor 1 would be obtained from attractor 
2 by the transformation 21 <-^ 23, 22 <-^ 24. 

For small values of a, the dominant oscillator of each of the attractors has a large amplitude, while the other two 
oscillators remain small. That means, the three attractors are located in phase space close to the subspace of the 
respective dominant oscillator together with the nonlinear trigger xj. Thus, they are clearly separate from each other. 
With increasing a, the attractors grow and begin to interact. First, attractor 3 vanishes at a ~ 0.052 in a boundary 
crisis, where it collides with the separatrix to the basin of attractor 1. Attractor 1 also loses stability at a ~ 0.06, 
such that only attractor 2 remains stable. Finally, at a = 0.11, a third crisis occurs, which leaves a large attractor 
that encompasses the dynamics of all three attractors. I n g eneral, the motion from one attractor to the other is much 
slower than the motion on the attractors (compare Fig. \U\ ). 

In Fig. O, the Lyapunov exponents as a function of a are shown. For a < 0.06, the GRS is on attractor 1. In 
the grey-shaded region a G [0.06,0.115], the GRS revolves on attractor 2. For a > 0.115, the Lyapunov exponents of 
the remaining hyperchaotic attractor can be seen. Here, the Lyapunov exponents grow approximately linear with a, 
similar to the case N = 5. 

C. Dynamics for arbitrary odd dimension A*' 

We have demonstrated how the dynamics of the GRS in the cases N = 5 and N = 7 can be understood within 
the same scheme of coexisting attractors. We expect that the dynamics of the GRS behaves in an analogous way for 
every odd dimension N. Thus, we expect, in general, ~ coexisting attractors for small values of a. On each of the 



attractors, one of the oscillators of the mode picture will be dominant. Each attractor will be located close to the 
hyperplane of the respective dominant oscillator together with the nonlinear trigger. The other oscillators will have 
a negligible amplitude. In the projection onto the dominant oscillator together with the nonlinear trigger xn, the 
dynamics will be Rossler-like. For higher values of a, the attractors will grow in phase space and interact with each 
other in several crises. Eventually, one hyperchaotic attractor will remain, which extends over all of the phase space 
that was occupied by the ^^^ coexisting attractors. 

For clarity, we shortly discuss the dynamics of the GRS in phase space for a hi ghe r number of variables, iV = 61. 



To this end, we use a pseudo space-time representation of the dynamics in Fig. 15 for three different values of a. 
The values of every second variable a;2n+i,?i = 0, . . . ,30, are depicted in grey scale as a function of time. In order 
to ensure that the dynamics has settled on the attractor, the GRS was propagated for a transient time 500000, 
which corresponds to about 10000 revolutions of the slowest oscillator for A'' = 61. A periodic orbit for a — 0.04 is 
visualized in Fig. ^q(a). In the mode picture, the mean amplitudes of the five dominant oscillators are: 1.30, 1.16, 
1.12, 0.40, 0.11, 0.05, 0.04. The amplitudes of the other oscillators are smaller than 0.01. Obviously, the dynamics 
is not dominated by a single oscillator, suggestin g th at the coexisting periodic orbits have undergone one or more 
merging crises before losing their stability. In Fig. g_5| (b), a comparatively low-dimensional chaotic orbit close to the 
above periodic one (Fig. pq(b)) is shown. It has been calculated for a = 0.07. The trace of the periodic orbit can 
still be seen in the image. Additionally, a considerable activity of high-frequency modes can be perceived as well. 
The largest Lyapunov exponents are: 0.00021, 0.00017, 0.00014, 0.00010, 0.00009, 0.00008, 0.00008, 0.00002, 0.00001, 
-0.00001, -0.00003, -0.00005. From this, the Lyapunov dimension is estimated to be £) = 20.8 with the help of 
the Kaplan- Yorke conjecture jlffl. In this state, the GRS does not explore all of its phase space and several other 
coexisting chaotic attractors can be found, which are not shown here. Figure |l^(c) shows a snapshot for a = 0.15, 
which lies in the hyperchaotic regime. Here, no coherent structure can be perceived. This state has 38 positive 
Lyapunov exponents and a Lyapunov dimension oi D = 60. In this state, the hyperchaotic attractor extends to all of 
the phase space and no coexisting attractor can be found. 

For large values of N, the GRS is similar to a time delay system of the form y{t) = f{y{t — t)) + g{y{t)), as 
investigated ,e.g., in [pT|-[l7[. This can clearly be seen in Fig. ^. The linear subsystem x acts as a delay-line, that 
merely transports signals back and forth between the positive fee dback process, represented by xi, and the nonlinear 



trigger xn- The delay time of the GRS is approximately N (see |III A| ) . The difference between the GRS and a time 
delay system is that the GRS - as a spatially discrete system ~ has got a nonlinear dispersion relation Lu{k). As every 
mode of the GRS is active for values of a in the hyperchaotic regime, the dispersion relation cannot be neglected in 
that case. Thus, the function / would have to be replaced by a functional in y(t — At) with Ai > N. Still, the GRS 
has got the same basic structure as a time delay system. We arrive at identical scaling properties of the Lyapunov 
exponents, the Lyapunov dimension, and the metric entropy for A'^ — > cxd, as will be shown in the following section. 

V. PROPERTIES OF THE GENERALIZED ROSSLER SYSTEM UNDER VARIATION OF SYSTEM 

SIZE N 

Up to now, we have investigated the dynamics of the GRS in phase space under variation of the control parameter 
a for the cases N — 5 and N — 7. In the present section, we investigate the GRS as a function of N for four fixed 
values of a: a — 0.3, a — 0.25, a = 0.15, and a — 0.07. The first three values of a are in the regime of large, phase-space 
filling hyperchaotic attractors, the fourth value of a lies in the regime of low-dimensional, coexisting attractors. N 
is varied from 3 to 61. Here, we mainly focus on the numerical calculation of Lyapunov exponents, from which the 
Lyapunov dimension Dx and the metric entropy H are estimated. The Lyapunov exponents have been calculated with 
the fixed initial condition (x, xn) = (1, 0, . . . , 0, 0). The other parameters are again chosen to be e = 0.1, b — 4, and 
d = 2. Interpreting the dimensionality N of the GRS as system size, we compare the results to the 'thermodynamic 
limes' of homogeneous, spatially extended systems, where the number of positive Lyapunov exponents, the Lyapunov 
dimension, and the metric entropy are reported to be proportional to N p8|-p2[. In this case, the distribution of 
Lyapunov exponents has been reported to approach a limit function /, i.e., Xi = f{i/V), if V is the system size 
p3[ . In the next step, we draw the attention to the GRS's similarity to time-delay systems, where the Lyapunov 
dimension has been reported to be proportional to the delay time, while the metric entropy approaches a limit value 
for increasing delay time UM- 

Figure ^ shows the number of positive Lyapunov exponents N^ , the Lyapunov dimension D\ (estimated via the 
Kaplan- Yorke relation flQ]) and the metric entropy H (estimated via the Pesin formula pj|) as a function of N for the 
four values of a. The number of positive Lyapunov exponents N~^ and the Lyapunov dimension D\ grow linearly with 
A'^. Remarkably, the Lyapunov dimension Dx is maximal up to its fractional part (i.e., Dx > N ~ 1) for almost all 
odd dimensions N, as long as the GRS is in the hyperchaotic regime (large enough a). The metric entropy approaches 



a limit value Hao for large N. Therefore, the Lyapunov exponents decrease proportional to 1/A^. Note that this 
observation is not so clear-cut for a = 0.07, probably due to the fact that the chosen initial condition lies within the 
basin of attraction of different attractors for different values of N. Therefore, we find that the limes N ^ oo oi the 
GRS is similar to the limes of an increasing delay time of time-delay systems and is not similar to the 'thermodynamic 
limes' of homogeneous, spatially extended systems. Recently, the constancy of the metric entropy of time-delay 
systems for increasing delay time has been attributed to the constancy of the number of 'localized nonlinearities' 
]l7[ . The same idea applies to the GRS. In Section II, we have distinguished between a linear, {N — l)-dimensional 
subsystem (xi, ..., a;(7v-i)) and a one-dimensional nonlinear subsystem (the trigger variable xn). The linear subsystem 
(x2, ...,X(jv-i)) (without the positive feedback process present in Xi) allows only for a linear, bi-directional transport 
of signals and is, therefore, not able to initiate any unstable behavior (see Fig. ^5|). In the case of the GRS, it is 
the one-dimensional 'localized nonlinearity' x^ |l7[ (together with xi), which allows for a nonlinear stretching and 
folding in phase space and is solely responsible for the positive metric entropy. As in the case of time-delay systems, 
we expect the value of the metric entropy H for increasing N to depend only on the number of localized nonlinearities, 
which is equal to one for the GRS, and the rate with which the information is processed, which is estimated via the 
correlation time Tc- Therefore, we hypothesize that the metric entropy scales like 

Hoc—, (10) 

Tc 

for N -^ oo. Since in the hyperchaotic case the correlation time is expected to be independent of N (see Sec. III. A), 
eq. (Kw reproduces the limit behavior of the metric entropy of the GRS. In the case of homogeneous, spatially extended 
systems, we expect the number of localized nonlinearities to increase proportional with the system size. This leads to 
a proportional increase of the metric entropy with the system size, as has been observed in several models. 

Thus, we conjecture that the high-dimensional chaotic motion observed in the GRS is fundamentally different 
compared to the spatio-temporal chaos observed in homogeneous, spatially extended systems. Here, the difference is 
expressed in terms of the number of localized nonlinearities a system possesses in phase space. In the case of the GRS, 
the number of the localized nonlinearities is independent of the control parameters including the system size. In the 
case of homogeneous and nonlinear spatially extended systems, the number of localized nonlinearities is expected to 
grow proportional to the system size. 

In Table fl we show the limit value Hoo for different values of a. Additionally, the value of Hoo normalized to the 
autocatalytic coefficient a is shown. Note that the value of a is an upper bound for the sum of the positive Lyapunov 
exponents and, with this, to the estimated metric entropy H . The case a = Hao corresponds to an unbounded, linear, 
{N — l)-dimensional system, (xi,a;2, ...^xn-i) in absence of the localized nonlinearity xn. 

From Table 1, we infer that the limit value of the metric entropy. Hoc, comes closer to its upper bound a for 
an increasing value of a, indicating that the hyperchaotic attractors increasingly exploit all the unstable directions 
available. 

In the following, we examine the distribution of Lyapunov exponents as a function of N . If one is to compare the 
Lyapunov exponents of different systems, one faces the problem that the absolute values of them depend on the chosen 
time scale. In systems with many different characteristic times, like the GRS, it is not obvious which time scale to 
use. Above, the time scale for the GRS under variation of N has been chosen such that the highest frequency of the 
linear subsystem ujmax is approximately equal to 2.0 (see eq. (pi)). In this case, we expect the correlation time Tc to 
be approximately constant for different values of N , while the time t^ a signal needs to traverse the linear subsystem 
(xi, ...,X(jv-i)) increases proportional to N. In the analysis to follow, we choose the time scale such that the rescaled 
mean positive Lyapunov exponent is equal to one. The rescaled Lyapunov exponents A^ are defined as 

N+ 
Xs ^ ^^A, (11) 

where X]A^ denotes the sum over the positive Lyapunov exponents, i.e., the metric entropy, and N^ denotes the 
number of positive Lyapunov exponents. Using the time scale according to eq. (pi]), the maximal frequency of the 
GRS scales like oJmax ~ 2^. Assuming N~^ x N and H — > Hao for high enough a, we find uJmax oc jj^. From 
this the correlation time is estimated to decrease with A^ in the hyperchaotic case like Tc oc -^^ . The same argument 
applies to the scaling of the velocity of signals yielding v ex -jj^ . This leads us to observe that the rescaling of the 
time according to eq. (^l|) makes the time r^ — N/v a signal needs to traverse the linear subsystem independent of 
the system size N for sufficiently large values of a. 

In Fig. n% the distribution of the rescaled Lyapunov exponents is shown for different values of A^ for a = 0.25. The 
Lyapunov exponents are sorted in descending order. The distributions converge to a limit distribution, X^^i = f{i/N), 
with increasing A^ . We observe qualitatively the same behavior for all values of a. Note that, although the existence 



of a limit distribution / has also been shown for homogeneous, spatially extended systems, the limes of an increasing 
system size iV for the rescaled Lyapunov exponents ( |ll|) is considerably different from the thermodynamic limes of 
spatially extended systems. In the latter case, the correlation time is approximately independent of the system size, 
while in the case of the GRS the correlation time decreases proportional to the system size N. 

In delayed dynamical systems with an expansive local dynamics (local in time), there has been reported the existence 
of "anomalous" Lyapunov exponents ||2^,|l^ which do not scale like 1/A^ for N -^ oo, but remain at a finite value. 
Considering the similarity of the GRS with a delayed system and the expansive term present in xi — axi — X2, one 
may expect to observe such an anomalous Lyapunov exponent for a > 0. However, due to the fact that the linear 
subsystem transports energy (in the form of squared amplitudes) away from xi, the local dynamics of the beginning 
of the linear subsystem is that of a damped oscillator for all values of a considered in this paper. For values of a being 
sufficiently large such that the local dynamics does indeed become expansive, the GRS is globally unstable, as the 
nonlinear trigger is no longer able to keep the amplitudes of the linear subsytem bounded. Accordingly, no sign of an 
anomalous Lyapunov exponent (i.e., a Lyapunov exponent that does not scale like 1/N) could be seen in any of the 
spectra calculated for the present work. 

Finally, we would like to discuss the shape of the limit distribution / of the rescaled Lyapunov exponents. We have 
observed above that the limit value of the metric entropy, if oc , approaches its upper bound a for an increasing value 
of a as a direct consequence of the dynamics increasingly exploiting all unstable directions of the unstable fixed point 
(x'^^^xJY ). We find that the same idea applies to the Lyapunov spectra: For an increasingly hyperchaotic attractor 
(increasing value of a), the Lyapunov spectra gain similarity with the real parts of the eigenvalues of the unstable 

fixed point {'x> \x)^ ). In Fig. 18, we show the Lyapunov spectra for three different values of a (fixed dimension 
N = 61) together with the real parts of the eigenvalues of the fixed point according to eq. (^. In all cases, the time 
scale has been chosen such that the trace of the matrix A equals one. 

VI. DISCUSSION AND CONCLUSION 

We have investigated the GRS as a model for high-dimensional chaos as the latter emerges out of low-dimensional 
chaos. One important feature of the GRS is that it consists of a linear subsystem with a variable number of degrees of 
freedom together with one nonlinear trigger. The linear subsystem can be solved analytically. Utilizing the eigenmodes 
of the linear subsystem, one can transform the GRS into a mode picture, consisting of harmonic oscillators that are 
coupled only via the nonlinear trigger. The mode picture reveals a structural symmetry of the GRS. With the aid 
of this structural symmetry, we interpret the dynamics of the cases N = 5 and N = 7 within a general scheme of 
coexisting attractors. For small values of a, there is, for each attractor, a specific projection into a three-dimensional 
subspace, where the dynamics of the GRS is Rossler-like. The attractors expand with increasing a and interact in 
several crises. This parameter regime of interacting attractors eventually leaves one large hyperchaotic attractor with 
many positive Lyapunov exponents. Even such kind of hyperchaotic dynamics can be made accessible to the human 
mind, which is used to envisage in three spatial dimensions, with the help of three-dimensional projections onto the 
oscillators of the mode picture. 

In the second part of this paper, we have investigated the Lyapunov exponents and related chaotic indicators of the 
GRS in the limit of large values of N, mainly in the hyperchaotic regime. The number of positive Lyapunov exponents 
and the Lyapunov dimension grow linearly with N. The Lyapunov dimension is maximal, Dx ~ N, independently of 
a, as long as one is in the hyperchaotic regime. The metric entropy converges to a limit value for increasing N. If 
the time is rescaled in such a way that the signal travelling time through the linear subsystem remains constant with 
increasing N, the distribution of the Lyapunov exponents approaches a limit function and the metric entropy grows 
linearly with N for N —>■ oo. We have argued that the hyperchaotic dynamics observed in the GRS is fundamentally 
different from spatio-temporal chaos. In this paper, we have expressed this difference in terms of the number of 
localized nonlinearities. 

In Section I, we have raised the question how dynamical systems develop from low-dimensional chaotic behavior 
to hyperchaotic states. In the GRS, one observes one specific path through the chaotic hierarchy, starting from a 
stable fixed point over chaos up to hyperchaos. The GRS exhibits a scenario of coexisting Rossler-like attractors 
that interact and eventually merge to form a hyperchaotic attractor. We would like to emphasize that the GRS only 
realizes one possible way, ending up with one special form of hyperchaos. We feel that there are many different forms 
of hyperchaos which possibly cannot be sufficiently characterized with the help of Lyapunov exponents. In the case 
of the GRS, the structure of the system has provided a helpful scheme to interpret the dynamics. We believe that 
the investigation of the topological structure of the fiow and, specifically, the attractor structure and the interaction 
of attractors, could be used to classify high-dimensional chaotic dynamics in general. 
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APPENDIX: NUMERICAL ALGORITHMS APPLIED 

The differential equations were integrated using a Runge-Kutta triple [g6|. This algorithm takes advantage of a 
sixth-order formula to propagate a system of ordinary differential equations. A fifth-order formula is used to estimate 
the integration error and, additionally, for each Runge-Kutta step t„ -^ in+i, the algorithm calculates a polynom 
that approximates the solution on the whole interval [t„,i„+i] up to an error of fifth order. The tolerance was set to 
10-10. 

The Poincar sections were obtained via parabolic interpolation in the vicinity of the intersection points. To calculate 
the bifurcation diagrams, we utilized a simple algorithm. For each value of a, we started with the state of the system 
for the last value of a, let the system adjust to the changed parameter value for a transient time equivalent to some 
thousand revolutions, and then recorded the intersection points. 

Lyapunov exponents were calculated using the algorithm described in [p7[-p9[. This algorithm tracks the time 
development of an orthonormal basis in the tangent space of phase space. At regular time intervals (AT = 50 in our 
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case), the vectors are reorthonormalized. The mean logarithmic growth rates of the moduH of the vectors arc the 
Lyapunov exponents. 

The delay equation of the time-delayed control was integrated by the same Runge-Kutta triple, using the spline 
polynoms that the Runge-Kutta triple outputs to record the continuous history of (x, xjv)- The control acts on all 
variables of the system. We choose a gain factor of 0.05 and limited the control signal to 25% of the modulus of the 
corresponding time derivative of the uncontrolled GRS. For the delay time r, we choose the eigenfrequency of the 
oscillator to be stabilized and, subsequently, adjusted r to minimize the mean control signal. 

Figure 1. Schematical representation of the structure of the GRS in (a) the Baier-Sahle picture, (b) the mode picture. Each 
circle represents one linear degree of freedom, the square represents a nonlinear degree of freedom. 

Figure 2. The real parts of the eigenvalues of the central fixed point {x.^^\xj^ ) close to the Hopf bifurcations that give rise 
to the coexisting attractors {N = 7) plotted as a function of parameter a. 

Figure 3. Bifurcation diagram for N — 5 under variation of a. Shown is the projection of the Poincare section onto 22 
(intersection at zi — 0, parameters; e = 0.1, & = 4, d — 2). 

Figure 4. Lyapunov spectrum for N = 5 under variation of a. The diagram shows the Lyapunov exponents A on attractor 1, 
except for the grey-shaded area, where attractor 1 does not exist (parameters: e = 0.1, fo = 4, d — 2). 

Figure 5. Period-1 orbit for N = 5 on attractor 1; 1:2 mode locking; projections onto (a) oscillator (21,22) and x^, (b) 
oscillator (23, 24) and x^ (parameters: a — 0.04, e = 0.1, 6 = 4, d = 2). 

Figure 6. Chaotic orbit for N = 5 on attractor 1. Projections onto (a) oscillator (21,22) and x^, (b) oscillator (23,24) and 
X5 (parameters: a = 0.085, e = 0.1, fe = 4, d = 2). 

Figure 7. Period-1 orbit for N=5 on attractor 2; 2:1 mode locking. Projections onto (a) oscillator (21, 22) and X5, (b) oscillator 
(23, 24) and x^ (parameters: a = 0.04, e = 0.1, b — 4, d = 2). 

Figure 8. Chaotic orbit for N = 5 on attractor 2; projections onto (a) oscillator (21,22) and x^, (b) oscillator (23,24) and X5 
(parameters: a = 0.1, £ = 0.1, fe = 4, d — 2). 



Figure 9. The transient onto the period-1 orbit of attractor 2 under the action of the time-delay control: (a) projection onto 
oscillator (21,22) and X5, (b) projection onto oscillator (23,24) and x^, (c) time development of the amplitude A of the control 
signal |x(f) — x(i — t)\ (parameters: N — 5, a — 0.12, e = 0.1, b = 4, d — 2; parameters of the control: r = 3.8557, control 
gain K = 0.05, control limit — 0.25). 

Figure 10. Time development of the mean amplitudes ri — \/ z\ -I- z\ and ri = y 2| -I- z\ of the two oscillators in the mode 
picture after the time-delay control has been switched ofi'. The trajectory remains in the basin of attraction of attractor 1 until 
t ~ 9500. At this point, it quickly moves into the basin of attraction of attractor 2 and never returns to attractor 1 afterwards 
(parameters: iV = 5, a = 0.12, e = 0.1, 6 = 4, d = 2). 

Figure 11. Hyperchaotic orbit for N — 5; projections onto (a) oscillator (21,22) and xs, (b) oscillator (23,24) and X5 
(parameters: a = 0.3, e = 0.1, 6 = 4, d = 2). 
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Figure 12. Bifurcation diagram for N — 7: attractor 1, attractor 2, and attractor 3 from bottom to top in different shadings. 
Tlie values of xe at the maxima of xe are shown (parameters: e = 0.1, 6 = 4, d = 2). 

Figure 13. Lyapunov spectrum for N — 7: depicted is up to a = 0.06 the development on attractor 1, from a — 0.06 to 
a — 0.11 that on attractor 2 (grey-shaded), above a = 0.11 that on the remaining large attractor. 

Figure 14. (a)-(c): Periodic orbit for N = 7, a = 0.035 on attractor 2 with 1:1:1 mode locking: projections onto (a) oscillator 
(21,22) and xr, (b) oscillator (23,24) and X7, (c) oscillator {z5,Zq) and xr. (d)-(f): Chaotic orbit for N = 7, a — 0.095 on 
attractor 2; projections onto (d) oscillator (21,22) and xr, (e) oscillator (23,24) and xj, (f) oscillator (25,23) and xr (other 
parameters: £ — 0.1, 6 = 4, d = 2). 

Figure 15. Time development of the GRS for A'^ = 61. The amplitudes of the variables X2n-i are shown in grey scale as a 
function of time for (a) a — 0.03, (b) a = 0.07, (c) a = 0.15 (other parameters: e = 0.1, b — 4, d = 2). 

Figure 16. (a) Number of positive Lyapunov exponents, (b) Lyapunov dimension, and (c) metric entropy as a function of A'^ 
for a = 0.3, a — 0.25, a = 0.15, a = 0.07 (other parameters: e = 0.1, b = 4, d = 2)). 

Figure 17. Distribution of rescaled Lyapunov exponents As = ^^^^ — A for different A*' (A'^ = 15,21,41,61). The Lyapunov 
exponents are sorted in descending order (parameters: a = 0.3, e = 0.1, 6 = 4, d = 2)). 

Figure 18. Comparison of the real parts of the eigenvalues of the linear subsystem a; and the Lyapunov exponents of the 
GRS in the case N = 61 for a = 0.07, a — 0.15, and a = 0.25. Here, all values have been divided by the respective values of a, 
in order to get comparable values for different a. (other parameters: e = 0.1, 6 = 4, d = 2). 



H^ 



Hex 



0.3 0.2 0.7 

0.25 0.16 0.6 

0.15 0.08 0.5 

0.07 0.02 0.3 

0.03 0.0 0.0 

Table L Limit value of the metric entropy H for different values of a. 
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